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Mechanism of generation of the emission bands 
in the dynamic spectrum of the Crab pulsar 
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ABSTRACT 

We show that the proportionately spaced emission bands in the dynamic spectrum of 
the Crab pulsar (Hankins T. H. & Eilek J. A., 2007, ApJ, 670, 693) fit the oscillations 
of the square of a Bessel function whose argument exceeds its order. This function has 
already been encountered in the analysis of the emission from a polarization current 
with a superluminal distribution pattern: a current whose distribution pattern rotates 
(with an angular frequency w) and oscillates (with a frequency fl > uj differing from 
an integral multiple of to) at the same time (Ardavan H., Ardavan A. & Singleton J., 
2003, J Opt Soc Am A, 20, 2137). Using the results of our earlier analysis, we find 
that the dependence on frequency of the spacing and width of the observed emission 
bands can be quantitatively accounted for by an appropriate choice of the value of 
the single free parameter Q/lo. In addition, the value of this parameter, thus implied 
by Hankins & Eilek's data, places the last peak in the amplitude of the oscillating 
Bessel function in question at a frequency {^ rt^/uj'^) that agrees with the position 
of the observed ultraviolet peak in the spectrum of the Crab pulsar. We also show 
how the suppression of the emission bands by the interference of the contributions 
from differring polarizations can account for the differences in the time and frequency 
signatures of the interpulse and the main pulse in the Crab pulsar. Finally, we put 
the emission bands in the context of the observed continuum spectrum of the Crab 
pulsar by fitting this broadband spectrum (over 16 orders of magnitude of frequency) 
with that generated by an electric current with a superluminally rotating distribution 
pattern. 
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1 INTRODUCTION 

Very soon after the discovery of pulsars, it was realized that 
the very stable periodicity of the mean profiles of their pulses 
could only result from a source that rotates, and which there- 
fore p osesses a rigidly rotating radiation distribution l|Goldl 
119681 ). In this paper, we show that this source rotation is 
not only responsible for the periodicity of the pulses, but 
also determines the detailed frequency dependence of the 
emitted radiation. By inferring the values of two adjustable 
parameters from observational data (values that are consis- 
tent with those of plasma frequency and electron cyclotron 
frequency in a conventional pulsar magnetosphere) , and by 
mildly restricting certain local properties of the source, we 
are able to account quantitatively for the emission spectrum 
of the Crab pulsar over 16 orders of magnitude of frequency. 
The rigid rotation of the overall distribution pattern of 



the radiation from a pulsar is described by an electromag- 
netic field whose distribution depends on the azimuthal an- 
gle ip only in the combination (ip—ut), where u; is the angular 
frequency of rotation of the pulsar and t is the time. As we 
show in Appendix A, Maxwell's equations demand that the 
charge and current densities that give rise to this radiation 
field should have the same time dependence. Therefore, the 
observed motion of the radiation pattern of pulsars can only 
arise from a source whose distribution pattern rotates rigidly, 
i.e. a source whose average density depends on ip only in the 
combination {(p — uit). Furthermore, if a plasma distribution 
has a rigidly rotating pattern in the emission region, then 
it must have a rigidly rotating pattern everywhere; we show 
in Appendix B that a solution of Maxwell's equations with 
the time dependence d/dt = —ud/dtp applies either to the 
entire volume of the magnetospheric plasma distribution or 
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to a region whose boundary is an expanding wave front that 
will eventually encompass the entire magnetosphere of the 
pulsar. 

Unless the plasma atmosphere surrounding the pulsar 
is restricted to an unrealistically small volume, it is there- 
fore an inevitable consequence of the observational data that 
the macroscopic distribution of electric current in the mag- 
netosphere has a rigidly rotating pattern whose linear speed 
exceeds the speed of light in vacuo c for r > c/ui, where r 
is the radial distance from the axis of rotation|j Although 
special relativity does not allow a charged particle with a 
non-zero inertial mass to move faster than c, there is no 
such restriction on the speed of propagation of the varia- 
tions in a macroscopic charge or current distribution. For 
instance, the distribution pattern of a polarization current 
with a propagation speed that exceeds c can be created by 
the coordinated motion of aggregates of particles that move 
slower than c (iBolotovskii fc Ginzbur3l972l : iGinzburgJI 19721 : 
iBolotovskii fc Bvkovl [l990|). Such a polarization-current 
density is on the same footing as the current density of free 
charges in the Ampere-Maxwell equation, so that its propa- 
gating distribution pattern radiates as would any other mov- 
ing s ource of electromagnetic fields (|Bolotovskii fc Bvkovl 
llQQOf ): indeed, such superluminal polarization currents have 
been demonstr ated to be efficient emitters of radiation in 
the laboratory (Bessarab ct al. 2004; Ardavan et al.ll2004bl : 
ISingleton eTal. 2004 ; Bolotovskii fc Scrov 2005] )! 

Once it is acknowledged that the electric cur- 
rent emitting the observed pulses from pulsars has a 
superluminally-rotating distribution pattern, results from 
the published li terature on the electrodynamics of super- 
luminal sources llArdavanl | l998l:lArdavan etahlbool l2004al . 
I2OO7I : ISchmidt et al.l I2OO7I : lArdavan et al.l l2008al lbh can be 
applied to pulsars. In this paper, we explain the recently- 
observed emi ssion bands in the dy namic spectrum of the 
C rab pulsar ([Hankins fc Eilekl 120071 ) using the calculations 
of lArdavan et al.l ( 2003 ) . Both the oscillations of the inten- 
sity and the frequency dependence of the spacing and width 
of these bands are described by a Bessel function with a 
single parameter that is already constrained by other ob- 
servational data on the spectrum of the Crab pulsar. This 
Bessel function is characteristic of the spectrum of the ra- 
diation by a superluminal polarization current whose distri- 
bution pattern rotates (with an angular frequency to) and 
oscillates (with a frequency fl > uj difi'ering from an integral 
multiple of ui) at the same time. It difi'ers from the Bessel 
function encountered in the analysis of synchrotron radia- 
tion only in the relative magnitudes of its argument and 
its order: while the Bessel function describing synchrotron 
radiation has an argument smaller than its order and so 
decays exponentially wi th increasing frequeii cy, the Bessel 
function encountered in lArdavan et al.l (|2003J), whose argu- 
ment exceeds its order, is an oscillatory function of frequency 
with an amplitude that decays only algebraically. It is this 
slower decay of the Bessel function in question that endows 
the emission from a superluminally-rotating source with a 



^ A corollary of this implication is that the magnetospheric 
plasma of a pulsar cannot be fully charge-separated, as is com- 
monly assumed in the wor ks based on the Goldreich- Julian model 
iJGoldreich fc JulianllT969h . 



broad spectrum. The physical mechanism underlying the 
broadba nd nature of this emission is focusing in the time 
domain iJArdavan et al.l |2003| . l2008al ) : contributions toward 
the intensity of the radiation made over an extended period 
of emission time are received during a significantly shorter 
period of observation time. 

This paper is organized as follows. Section [2] de- 
scribes how the most radiatively efficient parts of a pul- 
sar magnetosphere are thin filaments within the super- 
lumin ally rotating part o f its current distribution pat- 
tern JArdavan et al.l |2007|). A knowledge of the morphol- 
ogy of these filaments l Ardavan et al.ll2007l ). and the ability 
of a source that travels faster than its own waves to make 
multiple contributi ons to the signal receiv ed by an observer 
at a given instant (| Ardavan et al.ll2004ah . are necessary to 
understand many of the traits of pulsar observations de- 
scribed later in the paper. Sectio n O summarizes our ear- 
lier analysis (|Ardavan et al.l [2003J) on the frequency spec- 
trum of the radiation from a rotating superluminal source. 
Using these results, we derive various features of the emis- 
sion bands in Section |4] and compare them to the frequency 
b ands seen in the interpu lses of the Crab pulsar [Figs. 6- 
8. iHankins fc Eilekl ((20071)]; using a single input parameter, 
related to the pulsar's plasma density, we reproduce the ob- 
servational bands and predict the final ultra-violet emission 
peak of the Crab pulsar. Section [S] discusses the suppression 
of the bands by interference; this is relevant to the much less 
frequency-banded microburs ts and nanoshots of the Crab's 
main pulses [Figs. 2 and 3, iHankins fc Eilekl (|2007l )]. Sec- 
tion |6] describes the continuum spectrum of the Crab pul- 
sar; by introducing a further input parameter related to the 
plasma dynamics we are able to account quantitatively for 
the whole emission spectrum over 16 orders of magnitude of 
frequency. Section [7] gives a short discussion and summary. 
The mathematical details of our arguments are presented in 
Appendices A, B and C. 



2 THE EMITTING REGION OF A ROTATING 
SUPERLUMINAL SOURCE: MULTIVALUED 
NATURE OF THE RETARDED TIME 

At large distances from the source, the radiation field of a 
superluminally rotating extended source at an observation 
point P is dominated by the emissions of its volume ele- 
ments that approach P along the radiation direction with 
the s peed of light and zero acce leration at the retarded 
time l| Ardavan et al.ll2004al . 120071 ). These elements consti- 
tute a filamentary part of the source whose radial and 
azimuthal widths become narrower (as Sr ~ Rp~^ and 
5tp ~ Rp^'^, respectively), the larger the distance Rp of 
the observer from the source, and whose length is of the 
order of the len gth scale Iz of the s ource parallel to the 
axis of rotation ( Ardavan et al.ll2007h . (Here, r, ip, and z 
are the cylindrical polar coordinates of the source points.) 
For an observation point P with spherical polar coordinates 
{Rp,ipp,9p) that is located in the far zone, this contribut- 
ing part of the source lies at r = esc 9p,ip = ipp + 3n/2, and 
is essentially a st raight line parallel to the rotation axis, as 
shown in Fig. [l] (|Ardavan et al.ll2007l ). (The dimensionless 
coordinate f stands for ruj/c, where to, the angular frequency 
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Figure 1. Schematic illustration of the light cylinder r = c/u), 
the filamentary part of the distribution pattern of the source that 
approaches the observation point with the speed of light and zero 
acceleration at the retarded time, the orbit of this filamentary 
source, and the subbeam formed by the bundle of cusps that 
eman ate from the constitu ent volume elements of this filament 
[after lArdavan et alj l|2003)]- The subbeam is diffractionless in the 
polar direction. The figure represents a snapshot corresponding 
to a fixed value of the observation time tp . 



of rotation of the distribution pattern of the source, is the 
same frequency as that with which the pulsar rotates.) 

Once a source travels faster than its emitted waves, 
it can make more than one retarde d contribution to the 
field observed at any given instant (JArdavan et al.l |2003| . 
l2004al . l2008al ). This multivaluedness of the retarded time 
means that the wave fronts emitted by each of the contribut- 
ing elements of the source possess an envelope, which in 
this case consists of a two-sheeted, tube-like surface whose 
sheets meet tangentially alon g a spiraling cusp curve; see 
Fig. [2] IjArdavan et aLll2004al ). For moderate superluminal 
speeds, the field inside the envelope receives contributions 
from three distinct values of the retarded time, while the 
field outside the envelo pe is infiuenced only by a single in- 
stant of emission time l|Schmidt et al.l 120071 ). Coherent su- 
perposition of the emitted waves on the envelope (where two 
of the contributing retarded times coalesce) and on its cusp 
(where all three of the contributing retarded times coalesce) 
results in not only a spatial but also a temporal focusing of 
the waves: the contributions from emission over an extended 
period of retarded time reach an observer who is located on 
the cusp during a significantly shorter period of observation 
time (JArdavan et al.ll2008al ). 

The field of each contributing volume element of the 
source is strongest, therefore, o n the cusp of t h e env elope 
of wave fronts that it emits [see lArdavan et al.l (|2007r ) and 
references therein]. The bundle of cusps generated by the 
collection of the contributing source elements (i.e. by the 
filamentary part of the source that approaches the observer 
with the speed of hght and zero acceleration) constitutes 
a radiation subbeam whose widths in the polar and az- 



imuthal directions are of the order of S6p '^ Rp~^ and 
5ipp ~ Rp~^, respectively IjArdavan et al.ll2007l ). The over- 
all radiation beam generated by the source consists of a 
(necessarily incoherenlQ) superposition of such subbeams. 
This beam's azimuthal width is the same as the azimuthal 
extent of the source, and its polar width, arccos(l/f<) ^ 
\0p — 7r/2| ^ arccos(l/f>), is determined by the radial 
extent 1 < f< ^ f ^ r> of the superluminal part of 
the source (JArdavan et al.ll2007^ . This will be important in 
Sees. |4] and [S] the ipp dependence of the radiation intensity 
within the overall beam (i.e. what is observed as the main 
pulse, interpulse, and other components of the mean pulse) 
thus refiects the distribution of the source density around 
the cylindrical surface r = esc 9p from which the main con- 
tributions to the field arise. 

Since the cusps represent the loci of points at which the 
emitted spherical waves interfere constructively [i.e. repre- 
sent wave packets that are c onstantly disp ersed and recon- 
structed out of other waves (|Ardavanlll998l )] . the subbeams 
generated by a superluminal source need not be subject to 
diffraction as are conventional radiation beams. Neverthe- 
less, they have a decreasing angular width (i.e. are non - 
diffracting) only in the polar direction IjArdavan et al.ll2007l ) . 
Their azimuthal width 3(pp decreases as Rp~^ with distance 
because they receive contributions from an azimuthal extent 
Sip of the source that likewise shrinks as Rp~^ . They would 
have had a constant azimuthal width had the azimuthal ex- 
tent of the contributing part of the source been independent 
of Rp. On the other hand, the solid angle occupied by the 
cusps has a thickness 5zp in the direction parallel to the 
rotation axis that remains of the order of the height l^ of 
the source distribution at all distances (see Fig. 1). Con- 
sequently, the polar width 59p of the particular subbeam 
that goes through the observation p oint decreases as Rp~ ^ , 
instead of being independent of Rp (JArdavan et al.ll2007l ). 

Because it has a constant linear width parallel to the 
rotation axis, an individual subbeam subtends an area of the 
order of Rp, rather than Rp^ . In order that the energy flux 
remain the same across all cross sections of the subbeam, 
therefore, it is essential that the Foynting vector associ- 
ated with this radiation correspondingly decay more slowly 
than that of a conventional, spherically decaying beam: as 
Rp~^ , rather than Rp~^ , within the bundle of cusps that 
emanate from the constituent volume elements of the source 
and extend into the far zone (lArdavan et al. 2004a). This 
result, which also follows from the superposition of the 
Lienard-Wiechert fields of the co nstituent volume elem ents 
of a rotating superluminal source (lArdavan et al.ll2007l). has 
been dem onstrated experimentally ( Ardavan et al.l l2004bl : 
ISingleton et al. 2004i ). 

The fact that the observationally-inferred dimensions 
of the plasma structures responsible f or the emission from 
pulsars are less than 1 metre in size (JHankins et al.l 120031 ) 



^ The superposition of the subbeams is necessarily incoherent be- 
cause the subbeams that are detected at two neighbouring points 
within the overall beam arise from two distinct filamentary parts 
of the source with essentially no common elements. The inco- 
herence of this superposition would ensure that, though the field 
amplitude within a subbeam, which narrows with distance, decays 
nonspherically, the field amplitude associated with the overall ra- 
diation beam, which occupies a constant solid angle, does not. 
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Figure 2. (a) Cross section of the Cerenkov-like envelope (bold 
curves) of the spherical Huygens wave fronts (fine circles) emit- 
ted by a volume element S within the distribution pattern of an 
exten ded, rotating s u perlum inal source of angular frequency lo 
[after lArdavan et alj l l2004al )]. The point 5 is on a circle of ra- 
dius r = 2.5c/aj, or f = ru/c = 2.5, i.e. has an instantaneous 
linear velocity ru) = 2.5 c. The cross section is with the plane of 
S's orbit; dashed circles designate the light cylinder rp = c/u 
ifp = 1) and the orbit of S. (b) Three-dimensional view of the 
light cylinder, the envelope of wave fronts emanating from S, and 
the cusp along which the two sheets of this envelope meet tan- 
gentially. The cusp is tangent to the light cylinder in the plane of 
the orbit and spirals outward into the far zone. 



reflects, in the present context, the narrowing (as Rp~ and 
Rp~^, respectively) of the radial and azimuthal dimensions 
of the filamentary part of the source that approaches the 
observer with the speed of light and zero acceleration at the 
retarded time. Not only do the nondiffracting subbeams that 
emanate from such filaments account for the nanostructure, 
and so the brightness temperature, of the giant pulses, but 
the nonspherical decay of the intensity of such subbeams (as 
Rp~^ instead of Rp~^) explains why their energy densities 
at their source appear to exceed the energy densities of both 
the plasma and the magnetic field at the suface of a neutron 
star when estimated on the basis of the inverse square law 
l|Soglanov et al.ll2004l ). 



FREQUENCY SPECTRUM OF THE 
RADIATION FROM A ROTATING 
SUPERLUMINAL SOURCE 



In thi s section we summarize the results of lArdavan et al.l 
((2003) relevant to pulsars; the original intent of that paper 
was to calculate the spectrum of the radiation from a generic 
superlu minal source that has been implemented in t he lab- 
oratory (| Ardavan et al.ll2004bl : ISingleton et al.ll2004l ). This 
source comprises a polarization-current density j = dP/dt 
for which 



Pr,v,z{r,'P,z,t) = 



with 



(p = if ~ cut, 



Sr,ip,z{r, z) cos(m(p) cos{Qt), 

—n < if) ^ TT, (1) 



(2) 



where Pr,ip,z are the components of the polarization P in 
a cylindrical coordinate system (r, ip, z) based on the axis 
of rotation, s(r, z) is an arbitrary vector that vanishes out- 
side a finite region of the (r, 2) space, and tti is a positive 
integer. For a fixed value of f, the azimuthal dependence 



of the polarization ([T]) along each circle of radius r within 
the source is the same as that of a sinusoidal wave train 
with the wavelength 2-Kr/m whose m cycles fit around the 
circumference of the circle smoothly. As time elapses, this 
wave train both propagates around each circle with the ve- 
locity rtu and oscillates in its amplitude with the frequency 
n. This is a generic source: one can construct any distri- 
bution with a uniformly rotating pattern, Pr,^,z{r,if>,z) or 
jr,ip,z{r,'^,z), by the superposition over m of terms of the 
form Sr,ip,z{r,z,m)cos{m(fi). In the following discussion, we 
assume that the modulation frequency fl is positive and dif- 
ferent from an exact integral multiple of the rotation fre- 
quency u) [for the significance of this incommensurablity re- 
quirement, see lArdavan et al.l l|2003n ]. 

Although formulated in terms of a polarization current 
(for which there is manifestly no restriction on the propa- 
gation speed of the variations in the current density), the 
results that follow hold true for any current distribution 
whose density depends on ip as ip ~ ujt in r > c/uj. The 
electromagnetic fields 



E = -VpA" 



-^dA/dti- 



B = VpxA, 



(3) 



that arise from such a source are given, in the absence of 
boundaries, by the following classical expression for the re- 
tarded four-potential: 

yl''(xp,ip) = c'^ I A^xAtf{yi,t)5(tp-t~R/c)/R 



M = 0,---,3. (4) 

Here, {xp,tp) = (rp,ipp,zp,tp) and (x,t) = {r,(p,z,t) are 
the space-time coordinates of the observation point and the 
source points, respectively, R stands for the magnitude of 
R = xp — X, and fi = 1,2,3 designate the spatial com- 
ponents . A and j, of A^ and j** in a Cartesian coordinate 



system i Jackson 19991). 

In lArdavan et al.l (|2003l ). we first calculated the 
Lienard-Wiechert field that arises from a circularly moving 
point source (representing a volume element of an extended 
source) with a superluminal speed rco > c, i.e. considered a 
generalization of the synchrotron radiation to the superlu- 
minal regime. We then evaluated the integral representing 
the retarded field of the extended source ^ by superpos- 
ing the fields generated by the constituent volume elements 
of this source, i.e. by using the generalization of the syn- 
chrotron field as the Green's fu n ction for the problem [see 
equation (19) of lArdavan et al.l (|2003l )]. For a source that 
travels faster than c, this Green's function has extended sin- 
gularities arising from the constructive interference of the 
emitted waves on the envelope of wave fronts and its cusp 
(see Fig. El). 

Correspondingly, the fields E and B are given, in the 
frequency domain, by radiation integrals with rapidly oscil- 
lating integrands whose phases are stationary on the loci 
of the coherently contributing source elements: those source 
elements that approach the observer, along the radiation di- 
rection, with the speed of light and zero acceleration at the 
retarded time. For a radiation angular frequency nui that ap- 
preciably exceeds the angular velocity lj|j these integrals can 

^ When discussing pulsars in the following sections, we shall be 
considering radiation frequencies in the GHz-ultaviolet range. 
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be evaluated by the method of stationary phase to obtain 
the following expression for the electric field of the emitted 
radiation outside the plane of the source's orbit: 



E = K<^Eo + 2^E„ exp{~inipp) 
in which 



En 



if-^ exp{-i[n{Rp + fyr) - {Q/iu)i^p + fyr)]} 



Q^ = (-l 



XQ^QrQ:, +{m 



,n^ -n}, 



M4 



M- 



2 



Qr ^ fy — f<, n < TTRp/{f> — r<) 



(5) 

(6) 
(7) 
(8) 



Qr^ {2nRp/nyexp{-m/i), n :^ ■KRp/{f>~r<y ,{9) 



and 



with 



[s^J„_n/„(n) + is,,J'„_n/^(n)]e|| 
+ [(s^cos6lp - s^smOp)J„_Q/^(n) 
-isr cos OpJ'„_Q/^ (n)] ex , 



dz exp{mz cos 9 p)sr 



(10) 



(11) 



[see equation s (15) , 



(23), 
Here, 



(46) and (66)-(68) of 
{f,z;fp,zp) stand for 



lArdavan et aD ()2003l )]. 
{ruj/c,zuj/c;rpu}/c,zpuj/c), {Rp,6p,(pp) are the spher- 
ical polar coordintes of the observation point P, 
/i± = {Q,/uj) ± m, f^ < 1 and f> > 1 are the upper 
and lower limits of the radial interval in which the source 
densities Sr,ip,z are non-zero, J and J' are the Anger function 
and the derivative of the Anger function with respect to its 
argument, respectively, and ey = e^ x n/je^ x nj (which is 
parallel to the plane of rotation) and ex = nxe|| comprise a 
pair of unit vectors normal to the radiation direction n (e^ 
is the base vector associated with the coordinate z). The 
symbol {m -^ —m, Q, -^ — fi} designates a term exactly like 
the one preceding it but in which m and fi are everywhere 
replaced by —m and —SI, respectively. 

The radiated power per harmonic per unit solid angle 
is therefore given by 

{27T)-'cRp^\E„f 



dP„/dnp 



B7^)-'(cVl^')csc'ep|Q.|'Q^'|Q,|^ (12) 



where dilp denotes the element sin 6pd6pdipp of solid an- 
gle in the space of observation points. The contribution of 
the term {m —> —m, Jl —> —Q} in equation ^ has been 
ignored here because the Anger functions for —0, < turn 
out to be exponentially smaller than those for positive f2 
(see Appendix C). 

There is no difference between an Anger function, 
Ji^(x), and a Bessel function of the first kind, Ji/(x)i when u 
is an integer. Even for a non-integral value of u, the dif- 
ference between these two functions vanishes, as x~^i il 



whereas the rotation frequency of the Crab pulsar is Lu/2Tr 
30 Hz. The condition nw S> l^ will therefore always hold. 



X 2> 1. In the regime n S> 1, where the radiation frequency 
ncj appreciably exceeds the rotation frequency, therefore, 
the Anger functions Jn-ii/Lj{n) and J'n-n/uj{n) in Equa- 
tion (|10p can be replaced by the Bessel functions Jn-n/uj{n) 
and J'„_n/t^(n), respectively. If the radiation frequency nu) 
appreciably exceeds also the modulation frequency Q, then 
these Bessel functions can in turn be approximated by the 
following Airy functions: 






" \nj \nj CO 



(13) 



(14) 



where Ai' stands for the derivative of the Airy function with 
respect to its argument (see Appendix C). 
Hence, the radiated power is given by 



dP„ 
dilp 



Stt^u 



■CSC 6'p|Q,.| Q^' 



2x2/3 



xAi 



[srBii -I- (s^ cos 6p — Sz sin 9p)eA^^ 
2\V3o" 






i(s^e|| -f Sr cos 6lpex) 



/2y/=' n 



(15) 



in the regime where the radiation frequency (nw) apprecialy 
exceeds both the rotation (w) and modulation (Jl) frequen- 
cies of the source distribution ([l| . 

Asymptotic values of the amplitudes of 
Ai[-(2/n)^''^(n/a;)] and Ai'[-(2/n)^''^(!^/a;)] are respec- 
tively given by (2/n)-i''i2(t7/^)-i/4 ^nd {2lnf^^'^{^luSyl\ 
when n <C (n/o^)^, and by the constants Ai(0) and Ai'(O) 
when n > (fi/tj)^ (|Abramowitz fc Stegunlll970l ). Moreover, 
the quantity Q^3, which appears in equation (|15p . is inde- 
pendent of n (~ 20, /ui) if n <C /i+, decays as n~^ if n ^ /i+, 
and is of the order of n if n ~ m ^ 0,/u [see equation ((Tjl]. 

Depending on the relative magnitudes of n, m and Q,/uj, 
therefore, the amplitude of the radiated power dP„/df2p de- 
cays with the harmonic number n as n~'^|sr,^,z(n)|^, where 
the exponent o can assume one of the values 1/2, 2/3, 5/2 or 
8/3, and Sr,^,z{n) denotes the frequency dependence of the 
Fourier components of the source densities Sr,ip,z. The way 
in which Sr,^,z{n) decay with n, for large n, is determined 
by the distribution of the source density s(r, z) with respect 
to z [see equation (|lip ]. and so is model-dependent |j How- 
ever, the spacings between the emission bands predicted by 
equation (|15|l turn out to be essentially independent of the 
rate of decay of the amplitudes of these bands (see below). 

Having assembled the relevant equations, we now com- 
pare their predictions with established observational data 
on the Crab pulsar. 



■* Note that the electric susceptibility of the magnetospheric 
plasma (contained in the factor s) does not influence these re- 
sults, because it depends on th e source frequencies fJ.± id, not on 
the radiation frequency nu! [see lArdavan et al.l ll2008a) ] 
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Figure 3. The power dP„/dSlp radiated per harmonic per unit 
solid angle (in arbitrary units), into the signal linearly polarized 
parallel to the plane of rotation, versus the harmonic number n 
(lower axis) for Q./uj = 1.9 X 10'*, ay = 1.48 X 10"^ and 6|| = 
1.03 X 10*. The upper axis gives observation frequency uldJI-k 
for the case of the Crab pulsar [aj/(27r) ~ 30 Hz]. The spacing, 
the width, and the number of oscillations in this figure are the 
same as those of the bands seen in the interpulse data from the 
Crab puls ar over the frequ e ncy i nterval 8 < !^ < 10.5 GHz [see 
Figs. 6-8 MHankins fc Eilekl ||2007|) ]. 



4 PREDICTED FREQUENCY BANDS: 

COMPARISON WITH THE SPECTRA OF 
THE CRAB INTERPULSES 

One of the most remarkable features of the observational 
data from the Crab pulsar is the presence of well-defined fre- 
quency bands within the midpulses (Hankins fc Eilek 2007) . 
We now show that these are a natural consequence of a ro- 
tating superluminal source and use their spacing to predict 
the high-frequency emission of the Crab pulsar. 

When the absolute value of their argument appreciably 
exceeds unity, the Airy functions in equation psp are rapidly 
oscillating functions of the harmonic number n with peaks 
that spread further apart with increasing n (Fig.|3]). The rate 
of increase of the spacing between the peaks in question is 
determined solely by the value of f2/tj. We now derive an 
explicit expression for the spacing An of the maxima of the 
radiated power dP„/dJlp as a function of n and compare 
the result w ith th e proportionately spaced emission bands 
observed bv lHankin s & Eilek (2007). 

Within the range n/w <C rz <C (il/cj)'^ of harmonic 
numbers, where they are oscillatory, the Airy functions in 
equation (|15|) can be further approximated by 

Ai[-(2/n)'/^(f7/c^)] ~-K-^l''(nl2fl^'^(SllLo)-^l^coi.'^ (16) 

and 

Ai'[-(2/ri)'/^(f7/Lj)] ^■K-~^l'^(nl2Y^I^'^(SllLofl^w^^ (17) 

where 

V' = |(2/n)i''2(n/a;)3/2 - 7r/4 (18) 

(see Appendix C). In this regime, therefore, equation (|15p 
reduces to 

dP„ c^csc^gp - 2^ 2 / 2a; \ 1/2 I ^t, - |2 ,.,,■, 



in which 






V\\ = Sr cos tp — i 


— s^ sm ip 


(20) 


and 






V± = (s^ cos 6p ~ 


Sz sinOp) cos ijj - 


20, 
i — Sr cos 0p sin f/" (21) 

ntu 



represent the contributions towards the radiated field E„ 
from differing polarizations. 

Let us first consider the contribution towards the radi- 
ated power from the component of the field parallel to the 
plane of rotation. It follows from equation (|20p that 



=(2711 ){ 



20 



1^11 1 



where 



1 [2n/{nc^)ll/29{ 

711 = 2 arctan ^^ 



I — i2 • 2 ^^ ^ 

\Sr\ sm 711 H IS 

nuj 



V? cos^lW 



^(V'-Tll)}, 



'4 



(22) 



(23) 



with 9 and * denoting the imaginary part and a complex 
conjugate, respectively. Two successive oscillations of jPyl^ 
in the neighbourhood of any given n arise, according to equa- 
tion (I22p . from a change in the value of the argument of 
cos^ {ip — 7||) by TT. Therefore, the spacing between two suc- 
cessive peaks of this function is locally given by 



An — TT TT— (w - 

' an 



-7||) 



(24) 



an equation that yields An versus n for all nlf| 

If the source densities Sr.^,z have a weak dependence on 
n, then the high-frequency (n ^ 1) limits of equations 
and (|23l) assume the simple forms 



and 

711 ~ 

where 



|sr|^cos^(V'-7||), 



/2ny/2 a{sr4} 

Knuo) ISrP 

[2Q./{nLo 



(25) 



11/2 



a|[ + 6||(n — no) 



9{SrS*} 



and 



d_ 
dn 



Q{Sr.S 






(26) 



(27) 



(28) 



denote the first two coefficients in the Taylor expansion of 
9{srS^}/|sr| about a local value no of n. 

In this case, equation (I24|l . together with equations (|18|l 
and (|25p - (|28p . yields explicit expressions both for the radi- 
ated power \V\\\^ and for the spacing An between its suc- 
cessive peaks. These expressions are plotted in Figs. |3] and |4] 



^ At any point such equations may be readily translated into 
the frequency units Au versus u employed in Fig. 10 of 



I Hankins &: Eilelj l l2007f) by multiplying both An and n by ij/2iT : 
30 Hz, the rotation frequency relevant for the Crab pulsar. 
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over the interval 2 x 10* ^ n ^ 3.5 x 10**. In the case of the 
Crab pulsar, where the rotational frequency u/{2n) is ap- 
proximately 30 Hz, the interval of harmonic numbers shown 
in Figs. [3] and |4] corresponds to a frequency [u — naj/(27r)] 
interval coveri ng 6 to 10.5 GHz, i . e. to t he interval of frequen- 
cies over which iHankins fc Eilekl (|2007|) observe the emission 
bands. 

We have chosen a value (1.9 ± 0.1) x 10* for the more 
fundamental free parameter il/uu that renders the slope 
dAn/dn predicted by equation (|24p the same as the slope 
0.058 ± 0.01 of the line pass ing through the observational 
data IjHankins fc Eilekl 120071 ) even when the phase 7|| is 
zero|j We have then fitted the An intercept of the resulting 
line to that of the data, without changing its slope, by fixing 
the values ay = (1.48±0.05)xl0"^ and by = (1.03±0.1)xl0* 
of the parameters that appear in the expression for 7|| . 

The Airy functions in equation (|15|l stop oscillat- 
ing once the value of their argument falls below unity 
JAbramowitz fc StegunllieTOD . We shall see in Sec. Ig] that 
the value (1.9 ± 0.1) x 10'* of the parameter Q/lo thus 
implied by the data of iHankins fc Eilekl (12007J) places the 
last peaks of these functions at a frequency [~ {Q,/ij)^uj/TT] 
that agrees well with the position u ~ 10*^ Hz of the ul- 
trav iolet peak in the observed sp ectrum of the Crab pul- 
sar (|Lvne fc Graham-Smithll2006l ). The widths of the emis- 
sion bands shown in Fig. [3] also agree with the widths of 
those that are observed: the number of bands falling within 
the interval 2.7 x 10** ^ n ^^ 3.5 x 10* in Fig. g] is the 
same as the number of bands that occur within the corre- 
sponding frequency inte rval (8.25-10.5 GHz) in Figs. 6-8 of 
IHankins fc Eilekl l|2007l ). 

The contribution towards the radiated power from the 
component of the field perpendicular to the plane of rotation 
can be written, in analogy with equations (|22|l and (|23|l . as 



iri 



:(27x){ - \{s^cost 



Sz svaUp 



l2 . 2 

sin 7x 



^i i I . o p p r I 

H \sr\ cos 9p COS 7_L + Sin — Sz sin 6'p 

nu) '- 

\Sr\ cos 6p] cos (lp 



-11-)], 



(29) 



where 



7x = 2 arctan 



[2n/(nt^)]^/^Q{(3y cosflj: 



.)s 



A 



Bp-s^ sinej,|^-[2r2/{ncj)l|srP cos^ 0p 



(30) 



These expressions imply a band spacing with essentially the 
same characteristics as those discussed above, differing from 
those in equations (|22|l and (|23|) only in that Sr and s^ 
exchange their roles with Stp cos 9p — Sz sin 6p and s,p cos 6p , 
respectively. 

It will therefore be seen that the salient parameter re- 
sponsible for the frequency bands in the Crab pulsar is the 
ratio n/uj = (1.9 ± 0.1) x 10*. Using Lj/2n ~ 30 Hz relevant 
for the Crab yields a modulation frequency Q/2tt ~ 570 kHz. 
A plausible cause for this modulation is plasma oscilla- 
tions within the pulsar's atmosphere. A plasma frequency of 
570 kHz implies a free-electron density A'^_ ~ 4 x 10^ cm"'*, 
a value that is consistent with the inferred properties of the 
atmospheres of neutron stars. 



^ Note that the non-integral value of fl/ui can be approximated 
by an integer everywhere except in the factor sin(7rf2/[i)) that 
appears in equation ((7]l. 
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Figure 4. The spacing An (right axis) of the emission bands 
shown in Fig. |3] versus their mean harmonic number n (lower 
axis). The solid line represents the function given in equation 
(HU for n/ui = 1.9 X 10-*, ay = 1.03 x 10'' and fey = 1.48 x lO'^. 
The upper axis shows observation frequecy u = nu!/2iT relevant 
for the Crab pulsar; the left axis sho ws Aw. Data poin t s for the 
Crab pulsar are taken from Fig. 10 of JHankins fc Eilekl ||2007|') . 



5 SUPPRESSION OF THE EMISSION BANDS 
BY INTERFERENCE: APPLICATION TO 
THE MAIN PULSES OF THE CRAB 

We now turn to the main pulses of the Crab pulsar, in 
whic h the frequ ency-banded structure is much less appar- 
ent JHankins fc Eilck 2007 ). If the components of the elec- 
tric current within the emitting region (the region in which 
the rotating distribution pattern of the current approaches 
the observer with the speed of light and zero acceleration) 
are such that V^ and V^ simultaneously contribute towards 
the value of the field, destructive interference of the contri- 
butions with differing polarizations could result in the sup- 
pression of the emission bands. 

This can be easily seen from equation (|19|) if we note 
that, in this case, the radiated power for n 3> 1 is propor- 
tional to 

■^^ oc Q^'^n''^^'^[cos^{tp - 7l|) + i^{n) cos^(V' - 7-l)], (31) 

where 



K(n) 



s^ cos f p 



; sin 6f 



(32) 



is a coefficient that, like 7|| and y±, could be expanded into 
a Taylor series about a local value no of n. 

Figure [5] shows the n dependence of the right-hand side 
of equation (|3H) for a case in which Q^p decays as n~^, 
n/uj = 1.9 x 10* (as in Figs. |3] and E)), ay = 1.03 x 10*, 
fey = 1.48 x lO"'^ (as in Fig. H]), 7x = 7|| + 7r/2, and k can 
be approximated by 0.9-|-0.15(n/no — 1) with no = 2 x 10*. 
The emission bands are thus replaced by small-amplitude 
modulations of the radiation intensity that gradually die 
out. The behaviour of the radiated power shown in Fig. [5] 
over the interval 2.7 x 10* ^ n Sj 3.5 x 10* is consistent 
with those of the short-lived microbursts that are observed 
over the corresponding frequency interval (8.25-10.5 GHz) 
within the main pulses of the Crab pulsar [see Figs. 2-4 of 
IHankins fc Eilekl (|2007l )]. 

The radiation described by equation (|15|l would be el- 
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Figure 5. The radiated power dP„/dQp per harmonic per unit 
solid angle (in arbitrary units) versus n from a region of the mag- 
netosphere in which the contributions from differing polarizations 
are comparable in magnitude. The dependence of the right-hand 
side of equation l|31|l on n is shown for Q^ ~ n^^ , 7± = 7|| +''r/2, 
K = 0.9 + 0.15(n/no — 1), and the same values of C/loi, an, by, 
and no as in Fig. |4] The top axis of this Figure gives frequency u 
relevant for observations of the Crab pulsar. 



liptically polarized if either s^ or Sr is dominant, and lin- 
early polarized if Sz is dominant. The high degree of lin- 
ear polarization of the interpulse in the frequency range 
8.25 - 10.5 GHz (JMoffett fc Hankind [19991 ') suggests, there- 
fore, that the interpulse arises from a region of the mag- 
netosphere in which the component of the electric current 
parallel to the rotation axis dominates, and so the radia- 
tion is polarized in the direction of ex , while the main pulse 
arises from a region in which more than one polarization 
contributes toward the observed field. 

This explains not only the difference between the 
frequency structures of the main pulse and the inter- 
pulse, but also the difference between their temporal struc- 
tures, i.e. the fact that the intensity of the interpulse 
has a broad and smooth distribution in time, while that 
of the mai n pulse is composed of randomly distributed 
nanoshots (JHankins fc Eilekll2007h . It would be possible to 
identify the individual nondiffracting subbeams described 
in Sec. [2] only in the case of a source whose length scale 
of spatial variations is comparable with RJ,^ (in the case, 
e.g. of a turbulent plasma with a superluminally rotating 
macroscopic distribution). The overall beam within which 
the nonspherically decaying radiation is detectable would 
then consist of an incoherent superposition of coherent, non- 
diffracting subbeams with widely differing amplitudes and 
phases, that would be observed as randomly distributed 
nanoshots. Otherwise, the smooth transition between the 
adjacent filamentary parts of the source that generate neigh- 
bouring subbeams would result in an overall beam that is 
likewise smoothly distributed. The filamentary sources sam- 
pled from an emitting region in which a single component 
of the electric current dominates would not be sufficiently 
different from one another in either amplitude or phase for 
their associated subbeams to be distinguishable. 

The observed difference between the dispersion mea- 
sures of the main pulse and the interpulse (JHankins fc Eilekl 
120071 ) further supports the fact that these pulses arise from 



distinct regions of the magnetosphere. As we have seen in 
Sec.[2j an observer who is located at the polar angle 9p sam- 
ples, in the course of a rotation, the distribution of the den- 
sity of the electric current around the cylinder f = esc 9p. 
Thus the profile of the mean pulse reflects the azimuthal 
distribution pattern of the source. 



6 THE BROADBAND CONTINUUM 
SPECTRUM OF THE RADIATION 

Since they only involve a limited range of frequencies, the 
features of the band emission we have discussed above do not 
depend on the spectral distributions of the Fourier trans- 
forms Sr.ip,z of the source densities Sr,^,z sensitively. To put 
the emission bands in the context of the continuum spec- 
trum of the Crab pulsar, however, we would need both a 
more realistic model of the emitting current distribution and 
a more explicit description of the dependence of the density 
of emitting current on frequency. This notwithstanding, the 
frequencies across which the observed spectrum of the Crab 
pulsar undergoes sharp changes can be identifled with the 
frequencies at which the Airy function and the coefficients 
Qip and Qr in the spectrum of the radiation described by 
equation (|15|) have their critical points. We have already 
seen that, for the value Q/u) = (1.9 ± 0.1) x 1 0** im plied 
by the observational data of iHankins fc Eilekl (|2007|), the 
position Q,'^ /{tvlu^) of the last peak of the Airy function in 
equation (|15[) coincides with the ultraviolet peak in the ob- 
served spectrum of the Crab pulsar (Lvnc & Graham-Smith| 
l2006l ) (Fig.©. 

The coefficient Q^ in equation (|15p is independent of n 
(~ 2Q,/lj) when n ^ /i+, decays as n~^ when n ^ fi+, and 
is of the order of n when n c^ m ^ Q/cu [see equation(l7ll]. 
Hence, if the frequency mu!/{2n) of the azimuthal variations 
of the source in the Crab pulsar [see equation ([l])] falls in 
the terahertz range, i.e. m ~ lO'^^, then the order of mag- 
nitude of Q^ would change from Q/lo ~ lO'* to n -^ 10^^ 
across the terahertz gap in its observed spectrum. The cor- 
responding enhancement of the radiated power by a factor 
of order (AQ^)^ ~ 10^* is in fact consistent with observa- 
tional data (see below and Fig. |6]). A plausible cause for a 
frequency mijj/2n ~ 3 THz would be cyclotron resonance of 
free electrons in a magnetic field of around 10® G. A mag- 
netic field of this order at, or close to, the light cylinder, is 
con sistent with the spin-down pr operties of the Crab pul- 
sar (JLvne fc Graham-Smithll2006l l. 

Moreover, it can be seen from Equations ([8]) and (O 
that the coefficient Qr changes from being independent of n 
to decaying as n~^ " when n increases past TTRp/{f^ — r<)^, 
and so the observation point falls within the Fresnel zone. 
Given that the Crab pulsar is at the distance Rp ~ 7.1 x 
10^^ cm and has a light cylinder with t he radius c/lo ~ 1.6 x 
10* cm (JLvne fc Graham-Smithll2006l l. this would account 
for the corresponding steepening of the observed spectrum 
of the Crab pulsar at n ~ 3x lO'^® (Fig. [6]) if the radial extent 
of the emitting plasma observed from Earth has the value 
fy—f^~ 6.8 X lO"'^, i.e. 6.8 percent of the light-cylinder 
radius. 

These generic features of the spectrum of the radia- 
tion from a superluminal source, which do not sensitively 
depend on the detailed distributions of the Fourier compo- 
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Figure 6. The spectral distribution log(dPi,/df2p), predicted 
by equation I I33I I. versus logu for u = nu;/(27r) ~ 30n Hz ap- 
propriate for the Crab pulsar. The points show observational 
data (where available) of th e spectrum of the Crab pulsar 
IjLvne fc Graham-Smithll2006h . In the model, the recovery of in- 
tensity at the ultraviolet peak at ~ 10^® Hz is caused by res- 
onant enhancement due to the azimuthal modulation frequency 
mLo/2iT ~ 3 THz. The change in gradient at ~ 10^* Hz reflects a 
transition across the boundary of the Fresnel zone (see text). 



nents Sr.ip,^ of the source densities, are adequately described 
by the following simplified version of equation (|15|l : 



dPn 

dQp 



/ON 2/3 



1/3 J7 



(33) 



where S{n) stands for the dependence of the source density 
on n. For Q/lo = 1.9 x 10* and a Qr{n) that steepens by 
the factor n~ ' across n ~ 3 x 10 , equation p3p yields 



the continuum spectrum shown in Fig. [6] if we assume that 
S{n) dominates over Q^ everywhere, except across n = lO^'^ 
where Q^ increases by the factor 10^'* (due to the above- 
mentioned resonance with mij/2n ~ 3 THz), and decays 
as n"^/^ in n < lO" and as n'^^^^ in n > lO". Note the 
close correspondence, in Fig. [S] of equation (|33p with the 
observational data on the spectrum of the Crab pulsar over 
16 or ders of magnitude of frequency IjLvne fc Graham-SmithI 
l2006l ). 



7 DISCUSSION AND SUMMARY 

We must stress the model-independent nature of the results 
we have reported here. The only global property of the mag- 
netospheric structure of the Crab pulsar we have invoked is 
its quasi-steady time dependence: that the cylindrical com- 
ponents jr,ip,z {r, ifi,z;t) of the density of the magnetospheric 
electric current depend on ip only in the combination ip — Lut. 
Not only does this property follow from the observational 
data (as rigorously shown in Appendices A and B), but it 
is one that has been widely used i n the published literature 
on pulsars [see, e.g.. lMestel et al.l (119761 )]. 

Calculating the retarded field that is generated by the 
Fourier component associated with the frequency n of such a 
quasi-steady current distribution, we have derived a broad- 
band radiation spectrum with oscillations whose peaks are 



proportionately spaced in, and decay algebraically with, fre- 
quency. By inferring the values of the parameters Q/2n ~ 
570 kHz and muj/2TT ~ 3 THz from the observational data 
and mildly restricting certain local properties of the source 
densities Sr,^,z that appear in the derived expression we have 
quantitatively accounted for the following features in the ob- 
served spectrum of the Crab pulsar: 

(i) the spacing of the emission bands, 

(ii) the frequency at which the extrapolated spacing be- 
tween these bands would reduce to zero, 

(iii) the widths of the emission bands (their number 
within a given frequency interval), 

(iv) the salient features of the continuum spectrum (its 
sharp rise across the terahertz gap, its ultraviolet peak, the 
change by 1 of the value of its spectral index at 10^* Hz), 
and 

(v) the differences in the polarizations, dispersion mea- 
sures, and the time and frequency signatures of the main 
pulse and the interpulse. 

Given that the derived band spacing An increases as 
n"^'^ with the harmonic number n [equations (|18[) and 
l|24p]. we predict that t he repetition of the observations of 
iHankins fc Eilekl (|2007|) at higher frequencies would result 
in a dependence of the band spacing on frequency that has 
a local slope steeper than 0.058. There should be 12 bands 
within the frequency interval 10-20 GHz, for example, whose 
spacings would increase with frequency as dAn/dn ~ 0.086. 

The two significant input parameters required to simu- 
late the Crab pulsar's spectrum are a modulation frequency, 
Q/2tt ~ 570 kHz, tentatively attributed to plasma oscilla- 
tions of electrons in the pulsar's atmosphere and mu!/2TT ~ 
3 THz, possibly corresponding to cyclotron resonance of 
electrons in a magnetic field of around lO'' G. The latter 
gives rise to the resonant increase in spectral weight in the 
ultraviolet (Fig. [S| . 

At this point, it is worth considering briefly the in- 
fluence of the ratio n/u) on pulsar spectra in general. 
Given that typical features of the pulsar's spectrum scale 
as (Q/lo)^, we might expect that pulsars with slower ro- 
tational frequencies (and therefore small values of ui) will 
have observed spectral intensities weighted towards higher 
frequencies. By contrast, so-called millisecond pulsars (with 
large u) might be expected to have emission concentrated 
at lower frequencies than the Crab. These predictions seem 
to be borne out both by Geminga {Lj/2n ~ 4.5 Hz), which 
has emission peaked in the ultra-violet to gamma-ray end of 
the spectrum, and by millisecond pulsars such as 19374-21 
{Lu/2n ~ 642 Hz) and 1957-f20 {uj/2-k ~ 622 Hz), both of 
which show no emis sion in the GHz range but str ong pulses 
at radiofrequencies iJLvne fc Graham-Smithll2006l ). 

Furthermore, once it is acknowledged that (as shown 
in Appendices A and B) the current emitting the observed 
pulses from pulsars has a superluminally-rotating distribu- 
tion pattern, the results reported in the published literature 
on the electrodynamics of superluminal sources (lArdavarJ 
1 1998"; Ardavan ct al. '20031. [2004al . l2007l : ISchmidt et al.ll2007l : 
[Ardavan ct al. 2008a, b) can be used to explain the ex- 
treme value s of the giant pulse s' brightness temperature 
(~ 10''° K) (JArdavan et"aDl2007t ) . temporal width (-^ 1 ns) 
and source dimension (~ 1 m) (JArdavan et al.ll2008al ). as 
well as the unique characteristics of the average pulses' 
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polarization (their occurrence as concurrent 'ortliogonal' 
modes with swinging position angles and with nearly 100 
per cent linear or circular polarization) ijSchmidt et al.ll2007l ) 
and spectra (the range —4 to —2/3 of their spectral indices 
and the brea dth of their bandwidt h, from radio waves to 
gamma-rays) (JArdavan et al.ll2008a ). 

We conclude by re-emphasizing that all these features 
of the observed radiation are consequences solely of the fact 
that the emitting current depends on the coordinate ip in 
the combination [ip — ujt) . This constraint may be inferred 
directly from the simplest observation, that pulsars have a 
rigidly-rotating radiation distribution, and Maxwell's equa- 
tions (see Appendices A and B); it leads to a current dis- 
tribution with a superluminally rotating pattern at a radius 
r > c/uj responsible for the unique features of pulsar emis- 
sion. The explicit form of the current distribution plays a 
role only in accounting for the pulse-to-pulse variations of 
the features discussed above (i.e. for the variable 'weather' 
in the pulsar magnetosphere, not its stable 'climate'). The 
plasma processes by means of which the magnetic field of 
the pulsar couples the rotational motion of the central neu- 
tron star to the observationally implied rigid rotation of the 
distribution pattern of the emitting current have no direct 
bearing on the results reported in this paper. The salient fea- 
tures of the observational data can be understood in terms of 
the superluminal emission mechanism without a knowledge 
of the magnetospheric structure of the pulsar. 
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APPENDIX A: RIGID ROTATION OF THE 
RADIATION PATTERN IMPLIES A RIGID 
ROTATION OF THE SOURCE'S 
DISTRIBUTION PATTERN 

As reflected in the highly stable periodicity of the mean 
profiles of the pulses detected on Earth, the pulsar radiation 
field has a rigidly rotating distribution pattern on average. 
That is to say, the cylindrical components of the received 
radiation fields E and B depend on the azimuthal angle ip 
only in the combination ip> — cut: 



Er,,p,z{r,(p,z;t) = Er,^.z{r,ip - ujt,z,t), 
Br,^,z{r, ip, z; t) — Br,^,z(r, p — ujt, z,t). 



(Al) 
(A2) 



where (r, ip, z) are the cylindrical polar coordinates based 
on the axis of rotation, and ui is the angular frequency of 
rotation of the observed radiation pattern. An equivalent 
statement is that the radiation fields E and B have a quasi- 
steady time dependence: 



d_ d_ 

dt dip 



Ej- 



0, 



and 



d d 



(A3) 



(A4) 



for the general solutions of these partial differential equa- 
tions are given by the expressions on the right-hand sides of 
equations (|A1|) and (|A2|) . 
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In the Lorenz gauge, the electromagnetic fields appear- 
ing in equation ([Sjl are given by a four-potential A'^ that 
satisfies the wave equation 



V^A^ 



1 yA^ 

C2 9t2 



47r , 

= — / 



(A5) 



where fi = 1,2,3 designate the Cartesian spatial compo- 
ne nts A and .i of the potential A'^ and the current density 
j^ (|jacksonlll999l ) . The retarded solution to this equation in 
unbounded space is given by equation Q, i.e. by 



yl''(xp,ip) = c ^ / d^a;dij''(x,t)G(x,t;xp,tp), 



where 



G{x,t:xp,tp) = 



S{tp - t - R/c) 
R 



(A6) 



(A7) 



is the corresponding Green's function (|jacksonlll999l ). 

Employing equation (|3]) to write out Br,f,,z in terms 
of the cylindrical components Ar,^,z of the vector poten- 
tial, we see that the potential also has the time depen- 
dence expressed in equations (|A1|| - (|A4[I . It follows from 
Br — r^^dAz/d(f> — dA^/dz, for instance, that 



L'Br — —— — i^A-z "7^^^^. 

r dip oz 



(A8) 



where jC stands for the differential operator appearing in 
equations (|A3|I and (IA4I) : 



dt dip' 



(A9) 



This and the corresponding equations for _B^ and B^ show 



depend on tp and t only in the combination ip — cot. 

On the other hand, the cylindrical components of the 
vector potential are related to the cylindrical components 
jr,ip,z of the current density via the following spatial part of 
equation (|A6|) : 



Az 



r Ar dp Az At G 



cos(</?p 
cos((^p 



p)jr + sin(v3p - 
p)j^ - sin(v3p 

Applying the operator 
d d 






Cp = 



dtp d(pp 



(AlO) 



(All) 



to both sides of equation (|A10|) and making use of the fact 
that G depends on ipp and ip in the combinations ipp — (p 
and tp — t [equation (|A7|l ]. and so CpG = —CG, we find 
that the resulting equation can be cast in the form 



-^rp 

VAz,\ 


= c M rArAipAzAtG 


X 


"cc 
cc 


s{ipp - p)Cjr + ava{ipp - p)Cj^ ' 
s{ipp - ip)Cj^ - sin(v9p - p)Cjr 



(A12) 



by means of integrations by parts with respect to ip and t. 
Hence, CArp,,f,p,zp vanishes if and only if £jV,^,z is zero. 



It follows from equations (|A9|) and (|A12|I . therefore, 
that a necessary and sufficient condition for the rigid ro- 
tation of the distribution pattern of the radiation field. 



CB 



rp,'Pp,zp 



= 0, is the corresponding rigid rotation 
^jr,^,z = of the distribution pattern of the source den- 
sity. 



APPENDIX B: RIGID ROTATION OF THE 
SOURCE'S DISTRIBUTION PATTERN 
EXTENDS BEYOND THE LIGHT CYLINDER 

Our purpose here is to show that a solution of Maxwell's 
equations that has the quasi-steady time dependence 

i.e. has a rotating distribution pattern, applies either to 
the entire magnetosphere or to an expanding region whose 
boundary propagates at the speed of light. By consider- 
ing the initial-boundary-value problem for the wave equa- 
tion (|A5[) with Cauchy data satisfying equation (|B1|) . we 
establish that such a solution cannot be smoothly matched 
to other types of solutions of Maxwell's equations across a 
boundary that is confined to a localized region of the mag- 
netosphere. 

The change of variables 



Ci = r, 6 



^3 = V~^t, (,4 = ip{r,ip,z,t), 



(B2) 



where tp is, for the moment, to be considered an arbitrary 
function, transforms the ^-component of the wave equation 
(IA5I1 into 



^ i 



dAz 



+ {\Vi>\ 



d^Az 



+ 



1 / 2 



1 
9M, 



+2V', 



d^Az 



+ 2 



9^4 



2 +2^/), 



,^i 



d^Az 

d^Az 

d^Az 



969^4 Vr2^^ ' c2'^V 969^4 
/ 2 , 1 , \ dAz 47r . 



(B3) 



where tpt = dtp/dt, tpr = dtji/dr, etc. The Jacobian of the 
above transformation is 



ate, 6, 6, ^4) 



dt}) dtp 
dt dp ' 



(B4) 



d{r,z,ip,t) 

which remains non-zero as long as tp is not a function that 
depends on p and t as in 93 — u)t. 

Now suppose that at a given time t the surface 
tp{r,p, z,t) — represents the boundary of a region within 
which the solution to equation (|B3|I satisfies the symmetry 
(|B1|) . i.e. within which Az is a function of ^1, ^2 and ^3 only. 
At all points on this boundary, we then have 



(B5) 



dAz 
d^4 


= 

e4=o 


and 


1 d A dAz\ d^Az 

Ci aci ^"^ a^i ) ' 96" 




/ 1 cj''\ d^Az 47r , 



= 0, 



(B6) 



C4=0 
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where HB6P is the wave equation under symmetry (|B1|I ; (|B5|) 
and (|B6|) hold at the boundary of the region in question by 
virtue of holding inside that region. 

Hence, any other type of solution of the wave equa- 
tion that would match the symmetric solution across ijj = Q 
smoothljiJ has to be sought by solving the following initial- 
value problem: given that the Cauchy data on the hyper- 
surface ^4 = are those expressed in HB5|) and (|B6|I . what 
is the solution to the hyperbolic partial differential equa- 
tion (|B3|) beyond ^4 — 0? Note that since dA^/dS^A, van- 
ishes at all points of the hypersurface ^4 = 0, its derivatives 
9^A,/aCiC4, d'^A,/di2^A and d'^A,/dS,-i£,A, which are in di- 
rections interior to this hypersurface, also vanish at ^4 = 0. 

Thus, equation (|B3|l in conjunction with the data (|B5[) 
and (IB6I) demands that 



ivvi' 



,^t 



d^A, 



9^4 



0. 



(B7) 



C4=0 



There are two ways in which this requirement could be met: 
either ^ = is a characteristic surface of the wave equation, 
and hence the first factor in (|B7I) vanishes, or at ?/; = 
the derivative of dA:,/dS,i in the direction normal to the 
boundary is also zero. In the first case, the boundary of the 
domain in which symmetry (IBlfl is satisfied, ^ = 0, will 
consist of an expanding wave front that propagates at speed 
c. In the second case, the solution outside 7/) = will also 
be symmetric, according to the Cauchy-Kowalewski theorem 
[cf.|Courant & Hilbcrt (1962)]. This is because the extension 
(by means of a Taylor series) of the data into an integral 
strip next to the boundary ip = will yield a solution that is 
again independent of ^4. Since the above argument applies 
also to the new boundary of the region thus extended, it 
follows that Az will also be be independent of ^4 outside the 
surface i/)(r, tp, z, t) = 0, i.e. will be a function of (r, ip~ujt, z) 
throughout the magnetosphere. 

The corresponding results for Ar and A^p may be de- 
rived in the same way. 






C = 



- n 

LO \ CO 



n\-i/3 



(Cl) 



we can write the leading terms in their asymptotic expan- 
sions for large fi as 

J^(m + Cm'^') ^ (2/M)'/'Ai(-2^/\), (C2) 

and 

J'(p + Cm'/') ^ -(2/M)'/='Ai'(-2^/\) (C3) 

[ see e quations (9.3.23), (9.3.27) of lAbramowitz fc StegunI 
|1970|)]. Taking the limits n ^ f2/a; of the coefficients and 
arguments of the above Airy functions, we obtain the ex- 
pressions in equations (I13|l and (|14|) . 

For (fJ/oj) ^ n ^ (Q./uj)'^, the argument of the result- 
ing Airy function and its derivative that appear in equations 
(|13l) and (|14l) is large, so that we can further use the asymp- 
totic approximations 

Ai(-x) - TT-^/^x"'/' cosdx^'/^ - ^/4), (C4) 

and 



Ai'(-X) 



-1/2 1/4 . /2 3/2 

: TT ' X sm(3X''^ 



V4), 



(C5) 



[see equations (1 0.4.60 ) and (10.4.62) of 

lAbramowitz fc StegunI () 197(3 )] to replace them by the 
trigonometric functions given in equations (I16|l and (|17|) . 

Note that if Q is replaced by —fl, the argument of the 
above Airy functions assume a positive value, and hence the 
resulting Anger functions ^n+n/Ljin) and J'n+n/t^('T-) decay 
exp onentially, rather than algebr aically, with increasing n 
[see lAbramowitz fc StegunI (|l970l) ]. 



APPENDIX C: ASYMPTOTIC EXPANSION OF 
THE SPECTRUM FOR HIGH RADIATION 
FREQUENCIES 

In this Appendix, we evaluate the leading term in the asymp- 
totic expansion of the Anger functions that appear in equa- 
tion © for a radiation frequency (no;) that appreciably ex- 
ceeds both the rotation (uj) and modulation (Q.) frequencies 
of the source, though not necessarily the frequency mij of 
its spatial oscillations. 

In this regime, the Anger functions Jn-Q/cj('^) and 
J'n-n/ujin) can be approximated by the Bessel functions 
Jn -n/uiM and J ' n-n/^ (n), respectively [see equation (41) 
of lArdavan et al.l (|2003l )]. Once we cast these Bessel func- 
tions into the canonical forms J/iifi + Ca*' ) ^nd J' fi{fJ. + 
CfJ-^ ) by introducing the new variables 



^ A discontinuity in the value of dAz /d^4 that represents a field 
component is not permitted here, for it would entail the introduc- 
tion of surface charges and currents with infinite densities within 
the magnetosphere. 



